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We analyze the perturbative and parametric stability of the QCD predictions for the azimuthal 
asymmetries in heavy quark leptoproduction. At leading order, the cos tp asymmetry vanishes 
whereas the cos2ip one is of leading twist and predicted to be about 15% at energies sufficiently 
above the production threshold. We calculate the NLO soft-gluon corrections to ^-dependent lep- 
toproduction to the next-to-leading logarithmic accuracy. The soft-gluon approximation provides 
a good description of the available exact NLO results at Q 2 < m 2 . Our analysis shows that, 
contrary to the production cross sections, the cos 2ip asymmetry is practically insensitive to soft 
radiation for Q 2 < m 2 at energies of the fixed target experiments. We conclude that the cos 2ip 
asymmetry is well defined in pQCD: it is stable both perturbatively and parametrically, and in- 
sensitive (in the case of bottom production) to nonperturbative contributions. Measurements of 
the azimuthal asymmetries would provide an excellent test of pQCD applicability to heavy flavor 
production. 
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I. INTRODUCTION 



In the framework of perturbative QCD, the basic spin-averaged characteristics of heavy flavor hadro-, photo- and 



electroproduction are known exactly up to the next-to-leading order (NLO). During the last ten years, these NLO 
results have been widely used for a phenomenological description of available data (for a review see [1]). At the 
C*") ■ same time, the key question remains open: How to test the applicability of QCD at fixed order to heavy quark 
production? The problem is twofold. On the one hand, the NLO corrections are large; they increase the leading 
■ order (LO) predictions for both charm and bottom production cross sections by approximately a factor of two. 
Oh, For this reason, one could expect that higher-order corrections, as well as nonperturbative contributions, can be 
O ' essential, especially for the c-quark case. On the other hand, it is very difficult to compare pQCD predictions 
for spin-averaged cross sections with experimental data directly, without additional assumptions, because of a 
high sensitivity of the theoretical calculations to standard uncertainties in the input QCD parameters. The total 
uncertainties associated with the unknown values of the heavy quark mass, m, the factorization and renormalization 
scales, up and /i/j, Aqcd and the parton distribution functions are so large that one can only estimate the order 
j_i ■ of magnitude of the pQCD predictions for total cross sections at fixed target energies [2,3]. 

In recent years, the role of higher-order corrections has been extensively investigated in the framework of the soft 
gluon resummation formalism. For a review see Ref. [4]. Soft gluon (or threshold) resummation is based on the 
factorization properties of the cross section near the partonic threshold and makes it possible to resum to all orders 
in a s the leading (Sudakov double) logarithms (LL) and the next-to-leading ones (NLL) [5-7]. Formally resummed 
cross sections are ill-defined due to the Landau pole contribution, and a few prescriptions have been proposed to 
avoid the renormalon ambiguities [8-10]. Unfortunately, numerical predictions for the heavy quark production cross 
sections can depend significantly on the choice of resummation prescription [11]. Another open question, also closely 
related to convergence of the perturbative series, is the role of subleading logarithms which are not, in principle, 
under control of the resummation procedure [11,12]. 

For this reason, it is of special interest to study those observables that are well-defined in pQCD. A nontrivial 
example of such an observable is proposed in [13,14], where the single spin asymmetry (SSA) in charm and bottom 
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production by linearly polarized photons, + N — > Q + X[Q], was calculated. 1 It was shown that, contrary 
to the production cross section, the single spin asymmetry in heavy flavor photoproduction is quantitatively well 
defined in pQCD: it is stable, both parametrically and perturbatively, and insensitive to nonperturbative corrections. 
Therefore, measurements of this asymmetry would provide an ideal test of pQCD. As was shown in Ref. [15], the 
SSA in open charm photoproduction can be measured with an accuracy of about ten percent in the approved 
E160/E161 experiments at SLAC [16] using the inclusive spectra of secondary (decay) leptons. 

In the present paper we continue the studies of perturbatively stable observables and calculate the radiative and 
nonperturbative corrections to the azimuthal asymmetry (AA) in heavy quark leptoproduction: 

1(1) + N(p) - 1(1 -q) + Q( PQ ) + X[Q]( Px ). (1.1) 

In the case of unpolarized initial states and neglecting the contribution of Z— boson, the cross section of the reaction 
(1.1) may be written as 

A; = wpM [1 + (1 - " )2] Q2) +2(1 -» «<*■«*> 

+2 (1 - y) a A (x, Q 2 ) cos 2<p + (2 - y)y/l-y a I (x, Q 2 ) cos tp}. (1.2) 



The kinematic variables are defined by 

S = (£ + p) 2 , Q 2 = -g 2 



Q 2 



2p-q' 
4m 2 

Q 2 = xyS, p=^~. (1.3) 



V 1 Q 2 „„.o 4m2 

p ■ V 



In Eq. (1.2), gt (<Jl) is the usual j*N cross section describing heavy quark production by a transverse (longitudi- 
nal) virtual photon. The third cross section, a a, comes about from interference between transverse states and is 
responsible for the SSA which occurs in real photoproduction using linearly polarized photons [13-15]. The fourth 
cross section, <tj, originates from interference between longitudinal and transverse components [17]. In the nucleon 
rest frame, the azimuth ip is the angle between the lepton scattering plane and the heavy quark production plane, 
defined by the exchanged photon and the detected quark Q (see Fig. 1). The covariant definition of ip is 



cosy = , y , smv?= v - — - — I n-e, (1.4) 

r „ = £ ^ pvqa e , n ^ = e^ a(3 p u q aPQ p. (1.5) 

In Eqs. (1.3) and (1.4), m and mjv are the masses of the heavy quark and the target, respectively. 

In leading order pQCD, the cosy dependence vanishes, af orn (x, Q 2 ) = 0. For this reason, in this paper we restrict 
ourselves to the azimuthal asymmetry, A(p, x, Q 2 ), associated with the cos2y distribution: 

M n2] = d 3 a;jv(y = 0) + d 3 o-,jy(y = tt) - 2d 3 <7 lN (<p = tt/2) = ea A (x 1 Q 2 ) 
{P,XM ' ff<T lN {<p = 0) + d*a lN (<p = n) + 2&<T lN {<p = n/2) a T (x, Q 2 ) + sa L (x, Q 2 ) ' [ '> 

where e = — ^ — ^--r and d 3 aiiy(tp) = -r—-^^—(p,x,Q 2 ,tp). Note that the asymmetry defined by Eq. (1.6) is 
1 + (1 - yY dxdt^dy 

simply related to the mean value of cos 2(p: 



1 The well-known examples are the shapes of differential cross sections of heavy flavor production which are sufHcientfy 
stable under radiative corrections. 
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FIG. 1. Definition of the azimuthal angle tp in the nucleon rest frame. 



2tt 

/ dtp cos 2<p 
o 



d 3 a m 



A(p,x,Q 2 ) = 2(cos2ip)(p,x,Q 2 ), 



(cos2tp)(p, x, Q 2 ) 



dxdQ 2 dp 



(1.7) 



fd<p 



d 3 a m 



(p,x,Q 2 ,ip) 



o 



dxdQ 2 dip 



In this paper, we calculate the NLO corrections to the cos2ip asymmetry to the next-to-leading logarithmic 
accuracy (so-called soft-gluon approximation). Also we analyze the nonperturbative contributions to the AA due 
to the gluon transverse motion in the target. Our main results can be formulated as follows: 

• The azimuthal asymmetry defined by Eq. (1.6) is of leading twist; at energies sufficiently above the production 
threshold, it is predicted to be about 15% for both charm and bottom quark production. 

• The soft-gluon approximation provides a good description of the available exact NLO results on leptoproduc- 
tion in the region of relatively low virtualities, Q 2 <m 2 ; when Q 2 ^> m 2 , the quality of the NLL approximation 
becomes worse. 

• Contrary to the production cross sections, the cos 2<p asymmetry in azimuthal distributions of both charm and 
bottom quark is practically insensitive to radiative corrections at Q 2 < m 2 . This implies that large soft-gluon 
contributions to the (^-dependent and ^-integrated cross sections cancel each other in Eq. (1.6) with a good 
accuracy. 

• pQCD predictions for the cos 2ip asymmetry are parametrically stable; to within few percent, they are insen- 
sitive to standard uncertainties in the QCD input parameters: pa, pf, ^-qcd and in the gluon distribution 



• Nonperturbative corrections to the 6-quark azimuthal asymmetry due to the gluon transverse motion in the 
target are negligible. Because of the smallness of the c-quark mass, the analogous corrections to A(p, x, Q 2 ) 
in the charm case are larger; they are of the order of 20% at Q 2 < to 2 . 

We conclude that, in contrast with the production cross sections, the cos2ip asymmetry in heavy quark lepto- 
production, A(p,x,Q 2 ), is an observable quantitatively well defined in pQCD: it is stable, both parametrically and 
perturbatively, and insensitive (in the case of bottom production) to nonperturbative corrections. Measurements of 
the AA in bottom production would provide an ideal test of the conventional parton model based on pQCD. 

Concerning the experimental aspects, azimuthal asymmetries in bottom production can, in principle, be measured 
at HERA using the angular distributions of secondary (decay) leptons [15]. AAs in charm leptoproduction can also 
be measured in the COMPASS and HERMES experiments. Due to the relatively low c-quark mass, data on the 
D-meson azimuthal distributions would make it possible to clarify the role of subleading twist contributions. 



function. 
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The paper is organized as follows. In Section II we analyze the LO and NLO parton level predictions for in- 
dependent leptoproduction of heavy flavor in the single-particle inclusive kinematics. We check the quality of the 
soft-gluon approach against available exact results and discuss the region of applicability of the NLL approximation. 
Hadron level predictions for A(p, x, Q 2 ) are given in Section III. We consider in detail the pQCD contributions and 
nonperturbative corrections to the cos 2^ asymmetry at the HERMES, SLAC, COMPASS and HERA energies. 

II. PARTONIC CROSS SECTIONS 
A. Born level predictions 

At leading order, (D(a em a s ), the only partonic subprocess which is responsible for heavy quark leptoproduction 
is the two-body photon-gluon fusion: 

l*(q)+g(k g )^Qip Q ) + Q( P ^. (2.1) 
The -y*g cross sections, <r]? orn (k = T,L,A,I), corresponding to the Born diagrams are [18,19]: 

*? OTn (p, *) = [(1 - xf + x 2 + p(l-x- H/2)] In i±| 

- [(l-xf -x(2-3x) + (l-x)p]p}, 
„ ire%a em a s r 1 + 8 1 

*Bom (/3j £ ) = ^ j _ pln _±l_ +2{1 _ (2 . 2) 

3 

<r 7 BOTn (p,i) = o. 



*Bom (/M) = ^J^ P 2 ^ - 25 - p/2)lnl±| - (1 -x)(l- 2x/p)p}, 



In Eqs. (2.2), eQ is the quark charge in units of electromagnetic coupling constant and we use the following definition 
of partonic kinematic variables: 

s=(q + k g f, 



2q ■ k g ' 



Am 2 , Am 2 

^V 1 "— " = 7TcF- (2 ' 3) 

Note that the cos ip dependence vanishes due to the Q <-> Q symmetry which, at leading order, requires invariance 
under tp — ► tp + it [20]. 

The hadron level cross sections, <Jk{x, Q 2 ) (k = T, L, A, I), have the form 

l 

r f Ax x\ Q 2 

a k (x,£) = J dzg{z,p F )a k [—,-\ 1 k g = zp, ^ = ^2' 

where g(z,pp) describes gluon density in a nucleon N evaluated at a factorization scale pf- The partonic cross 
sections, a k , are functions of p and x defined by Eq. (2.3). 

B. Soft-gluon corrections at NLO 

To take into account the NLO contributions, one needs to calculate the virtual <D(a em a 2 ) corrections to the Born 
process (2.1) and the real gluon emission: 
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7*(«Z) + 9(h) -» <3(pq) + Q(Pq) + fl(P fl )- 
The partonic invariants describing the single-particle inclusive (1PI) kinematics are 



(2.5) 



s' = 2q-k g = s + Q 2 = zS', 

S4 = s' + ti + Ml, 



h = (kg-PQY 

ui = (q-pq) 2 



m 

2 



m 



= zTi, 
Ui, 



(2.6) 



where s 4 measures the inelasticity of the reaction (2.5). The corresponding 1PI hadron level variables describing 
the reaction (1.1) are 



S' = 2q ■ p = S 
S 4 = S , + Ti 



Ui, 



Ti = (p-Pq) 

Ui = [q-PQ) 



(2.7) 



We neglect the photon-(anti)quark fusion subprocesses. This is justified as their contributions vanish at LO and 
are small at NLO [21]. 

The exact NLO calculations of the unpolarized heavy quark production in 717 [22,23], 7*5 [21], and gg [24,25] 
collisions show that, near the partonic threshold, a strong logarithmic enhancement of the cross sections takes place 
in the collinear, p g _T — > 0, and soft, p g — > 0, limits. This threshold (or soft-gluon) enhancement has universal 
nature in the perturbation theory and originates from incomplete cancellation of the soft and collinear singularities 
between the loop and the bremsstrahlung contributions. Large leading and next-to-leading threshold logarithms can 
be resummed to all orders of perturbative expansion using the appropriate evolution equations [5-7]. The analytic 
results for the resummed cross sections are ill-defined due to the Landau pole in the coupling strength a s . However, 
if one considers the obtained expressions as generating functionals of the perturbative theory and re-expands them 
at fixed order in a s , no divergences associated with the Landau pole are encountered. 

Soft-gluon resummation for the photon-gluon fusion has been performed in Ref. [26] and checked in Refs. [27,14]. 
To NLL accuracy, the perturbative expansion for the partonic cross sections, d 2 tJfc/dtidui (k = T,L,A,I), can be 
written in a factorized form as 



d 2 <Jfc 



(«', ti, ui) = Bf orn (s', ti,m)\ S(s' + h+ Ul ) + J2 



dtidui 

with the Born level distributions B® oln given by 
B^ orn {s',t 1 ,u 1 )=ne 2 Q a em a 
B£ or, V,ti,ui) =ire 2 Q a em a 
B5 orn (s',ii,ui) =ne 2 Q a em a 
Bf orn (s',£i,«i) =ire 2 Q a em a 



71=1 



Oi s C A 



K^\s',t u ux) 



t\ Ml 

— + -r + 4 

s' 



2 / 



m 2 s' 

txUi 



s'(m 2 - Q 2 /2) 
iiiti 



+ 



till! 



4V0 2 



2 / 

m s 

t\U\S 



2 / 
m s 



— + «! 

tmi s' 



m 



1/2 



h - Mi 



2Cf 



2m 2 s' 



(2.8) 



(2.9) 



Note that the functions (s' , ti, Ui) in Eq. (2.8) originate from the collinear and soft limits. Since the azimuthal 
angle ip is the same for both 7*5 and QQ ccnter-of-mass systems in these limits, the functions K( n \s', t\, u\) are 
also the same for all k = T,L, A, I. At NLO, the soft-gluon corrections to NLL accuracy in the MS scheme are 



A^V.ii.ui) 



In (s 4 / 



in 



s 4 



1 + ln 



I Mi 



- 1 



2C F 
C a 



(l + ReL /3 )+ln 



-Ml 



In 



(2.10) 
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where we use n = up = Hr. In Eq. (2.10), Ca — N c and Cf — (N 2 — 1)/(2N C ), where N c is number of colors, 
while Lp = (1 — 2m 2 /s){ln[(l — + /3)]+i7r}. The single-particle inclusive "plus" distributions are defined by 



V (s A /m 2 )' 


= lim \ 


s 4 





s{^»<»->+t> +, (^4 



For any sufficiently regular test function /i(s 4 ), Eq. (2.11) gives 



°4 

/ 



ds4 /l(s 4 ) 



V (s 4 /m 2 )" 


s 4 


s 4 





r , N , N1 ln Z (s 4 /m 2 ) 
ds 4 [/i(s 4 ) - /i(0)] — + 



s 4 



i + 1 



/i(0) In' 



Vm 2 



(2.11) 



(2.12) 



In Eq. (2.10) , we have preserved the NLL terms for the scale-dependent logarithms too. We have checked that 
the results (2.9) and (2.10) agree to NLL accuracy with the exact 0(a em a 2 ) calculations of the photon-gluon cross 
sections <tt and &l given in Ref. [21]. 

To perform a numerical investigation of the results (2.9) and (2.10), it is convenient to introduce for the fully 
inclusive (integrated over t\ and ui) cross sections, &k (k = T, L, A, I), the dimensionless coefficient functions cj 



(n,l) 



2 /. 2\ oo n / 2 \ 

n=0 1=0 ^ ' 



where the variable r\ measures the distance to the partonic threshold: 

s 



V = 



Am 2 



-1, 



m z 



(2.13) 



(2.14) 



Concerning the NLO scale-independent coefficient functions, only c^' ^ and c^ 1 ' ^ arc known exactly [21,28]. As 
to the ^-dependent coefficients, they can by calculated explicitly using the renormalization group equation: 

1 

da k (s',Q 2 ,n 2 ) 



dln/j 2 



= - J Az a k (zs' ,Q 2 , ^ 2 )P gg (z), 



(2.15) 



where z m j n = (4m 2 + Q 2 )/s' 7 &k(s' ,Q 2 , /z) are the cross sections resummed to all orders in a s and P gg (z) is the 
corresponding (resummed) Altarelli-Parisi gluon-gluon splitting function. Expanding Eq. (2.15) in a s , one can find 
[26,14] 

1 

4 M V,o = ^2 J dz -*)- (2-i6) 

2 m in 

where &2 = (HCU — 2n/)/12 is the first coefficient of the /3(a s )-function expansion and n/ is the number of active 
quark flavors. The one- loop gluon splitting function is [29]: 

1 



+ 



+ z(l -z))6(l-z-e) + 5(1 - z) lne \ C A + b 2 6(l - z). 



(2.17) 



With Eq. (2.16) in hand, we are able to check the quality of the NLL approximation against exact answers. In 
Figs. 2 and 3 we plot the functions c^' l \rj,£) and c^' l \rj,^) (n,l = 0, 1) at £ = 10~ 2 and £ = 3.16, respectively 
Predictions of the NLL approximation (2.10) are given by dotted curves. The available exact results are given by 
solid lines. One can see a reasonable agreement up to energies 77 w 2. As to the Q 2 -dependence, we have found 
that the soft-gluon approach reproduces satisfactorily the exact results at £ < 1. At high values of £, Q 2 3> m 2 , the 
quality of the NLL approximation becomes worse. 2 



2 Our analysis shows that the same situation takes also place for the energy and Q 2 behavior of the functions c^ 1 ''' (77, £), 
I = 0,1. We do not give corresponding plots for c^ 1 ''' (77, £) because the contribution of the longitudinal cross section to the 
cos 2<p asymmetry is small numerically. 



6 



0.001 0.01 0.1 1 10 100 1000 



1 








0.8 


\= io- 2 




-(0.0) 
C T 


0.6 








0.4 




(1,0) 
<-T 




0.2 


















0.001 0.01 0.1 1 10 100 1000 



0.001 0.01 0.1 1 10 100 1000 



0.2 


\= io- 2 


/ \ r 10 ' 01 
/ \ 


0.15 






0.1 






0.05 





' ""• r (1 ' 0) \ 
'. A \ 



0.001 0.01 0.1 1 10 100 1000 

T| = s I A m 2 - 1 



FIG. 2. c T ' (r?, £) and c A ' (rj, £) coefficient functions at £ 
the NLL approximation (dotted lines). 
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the NLL approximation (dotted lines). 
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III. HADRON LEVEL RESULTS 



A. pQCD predictions 

Let us now analyze the impact of the approximate NLO perturbative corrections on the AA at hadron level. We 
will consider the parameters A(Q 2 ), A(x) and A(y), 

2 _ 2 \^ cos2ifi wk M2M _ 2 ^ cos2 ^S (w) 



dQ 2 d(p yr, ~* ,r; o ~ r dxdip 



o r a o d aiN I \ 

2 J dycos2y) {p,y,<P) 

A ^ = 2^ 32 ' C 3 - 1 ) 

o dyd^s 

which describe the dependence of the cos2<^ asymmetry on Q 2 , Bjorkcn x and y, respectively. Unless otherwise 
stated, the CTEQ5M [30] parametrization of the gluon distribution function is used. The default values of the 
charm and bottom mass are m c = 1.5 GeV and m b = 4.75 GeV, A 3 = 260 MeV and A 4 = 200 MeV. For the 
factorization scale we use fiF = \/ m \ + 2 /4 in the case of bottom production and p,p = 2yjm 2 + Q 2 /4 in the 
charm case [1]. 

Our results for the Q 2 - and x-distributions of the AA in charm leptoproduction at several values of initial energy 
are presented in Fig. 4 and Fig. 5, respectively. The LO and NLO predictions are given by solid and dotted 
lines, correspondingly. The lines with label "1" correspond to p\ = 0.2, "2" — > p 2 = 0.1, "3" — > p3 = 0.05 and 
"4" -» p 4 = 0.025, where p = Am 2 /S. So, in the charm case, we have: E x = 24 GeV, E 2 = 47 GeV, E 3 = 95 GeV 
and E 4 = 190 GeV, where E is the lepton energy in the lab (nuclcon rest) frame: E = (S — m^-)/(2mjv)- The 
y-distribution of the asymmetry at the COMPASS energies is given in Fig. 6. In Fig. 7 we plot the A(Q 2 ), A(x) and 
A(y) distributions of the asymmetry in bottom production at the same values of pi = {0.2,0.1,0.05,0.025} which 
correspond to the following set of initial energies: Ei = {240, 480, 960, 1920} in units of GeV. 

Our calculations given in Figs. 4—7 represent the central result of this paper. One can see from Figs. 4 and 7 that 
soft-gluon corrections to A{Q 2 ) are about few percent at not large Q 2 < m 2 . At fixed values of x, the kinematical 
restriction Q 2 < m 2 leads to x < m 2 /S. For this reason, radiative correction to A(x) are small in the region of 
x < p (see Figs. 5 and 7). For comparison, we plot in Fig. 8 the so-called if-factors for (^-integrated cross sections: 

d(T NLO / d<r LO d<r NLO / d<7 LO 

K(Q ) = dQ'i / d g" and K{x) = ^ / d l ™ . One can see that large soft-gluon corrections to the production 
cross sections practically (to within few percent) do not affect the Born predictions for the cos2<p asymmetry at 
Q 2 < m 2 and x < p. 

At fixed values of y, the allowed region of Q 2 is mfy 2 /(l — y) < Q 2 < yS — 4m 2 , where m; is the initial lepton 
mass. Since production cross sections rapidly vanish with growth of Q 2 , practically whole contribution to A(y) 
originates from the low-Q 2 region. For this reason, radiative corrections to A(y) are negligible practically in the 
whole region of y (see Figs. 6 and 7). 

Let us now discuss the region of applicability of the adopted soft-gluon approximation. As noted in previous 
Section, soft radiation reproduces satisfactorily the existing exact NLO results when £ is not large, Q 2 < to 2 . At 
large Q 2 ^> to 2 , the hard (p s ,t 0) contributions becomes sizeable and the quality of the NLL approximation 
becomes worse for both (^-dependent and (^-independent cross sections. Moreover, we have observed that soft- 
gluon approximation overestimates the exact results for c^ ,0 \rj, £) and Cy 1 ' 1 ^ (77, £) at £ >• 1 and, simultaneously, 
underestimates the corresponding ones for c^' 1 -* (77, £). For this reason, in the high-Q 2 region, the full NLO corrections 
to the cos 2ip asymmetry may be smaller than the soft-gluon ones. 

As to the energy dependence, one can see from Figs. 2 and 3 that soft radiation describes very well the exact 
NLO results on ip- independent photon-gluon fusion at partonic energies up to r\ w 2. Since the gluon distribution 
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function supports just the threshold region, the soft-gluon contribution dominates the photon-hadron cross sections 
approximately up to S/Am 2 ~ 10 (and, correspondingly, up to p <~ 0.1 for <jin)- Using the exact results for the 7*5 
cross sections [28] , we have verified that the contribution originating from the region 77 > 2 makes only few percent 
from the NLO hadron-level predictions for a T (S, Q 2 ) and a L (S, Q 2 ) at S/Am 2 < 10 (and Q 2 < m 2 ). Results of Ref. 
[26] on soft-gluon corrections to F| harm (x,£) confirm our conclusion. 
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FIG. 4. Azimuthal asymmetry, A(Q 2 ), in c-quark production for several values of initial energy: Ei = 24 GeV, E2 = 47 
GeV, E3 = 95 GeV and E4 = 190 GeV. Plotted are the results at LO (solid curve), at NLO to NLL accuracy (dotted curve) 
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FIG. 5. Azimuthal asymmetry, A(x), in c-quark production. The notation is the same as in Fig. 4. 
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Azimuthal asymmetry, A(y), in c-quark production. The notation is the same as in Fig. 4. 
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FIG. 7. Asymmetry parameters A(Q 2 ) (left panel), A(x) (right panel) and A(y) (down panel) in fc-quark production for 
several values of initial energy: E ± = 240 GeV, E 2 = 480 GeV, E s = 960 GeV and E 4 = 1920 GeV. Plotted are the results 
at LO (solid curve), at NLO to NLL accuracy (dotted curve) and at LO with kr smearing (dashed curve). 
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FIG. 8. AT-factors in charm production at NLL level: K{Q 2 ) = / ^« (left panel) and K(x) = —± 

(right panel) . The notation is the same as in Fig. 4. 



Presently, exact NLO calculations of the (^-dependent cross section of heavy flavor production are not completed. 
However we can be sure that, at energies not so far from the production threshold, the soft radiation is the dominant 
perturbative mechanism in the case of (ja(S, Q 2 ) too. First, LO predictions for the cos 2</?-dependent cross section 
are large and the Sudakov logarithms have universal, ^-independent structure. For this reason, <ja(S,Q 2 ) has also 
a strong threshold enhancement. Second, our analysis of the exact scale-dependent cross section c^' 1 ^ , £) given 
in Figs. 2 and 3 confirms with a good accuracy the dominance of the soft-gluon contribution. These facts argue 
that hard and virtual corrections to the cos 2</?-dependent cross section cannot affect significantly the soft-gluon 
predictions for the azimuthal asymmetry at low Q 2 in the energy region up to S/4m 2 <~ 10. 
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FIG. 9. QCD LO predictions for A(£) in charm and bottom production at p — 1CP 3 (left panel) and p = 1CP 4 (right 
panel) with and without the inclusion of the k T smearing effect. 



Let us briefly discuss the origin of perturbative stability of the cos 2^3 asymmetry. Note that the mere ip- 

indcpcndcnt structure of the Sudakov logarithms cannot explain our results since perturbative stability does not 

c (i>o) c (o,o) 

take place at the parton level. In fact, the ratios -7^57(77, £) and -^oj(r], £) differ essentially from each other even 

at r),£ < 1. This is due to the fact that the physical soft-gluon corrections (2.8) are determined by a convolution 
of the Born cross sections with the Sudakov logarithms which, apart from factorized <5(s4)-terms, contain also non- 
factorizable ones (see Eq. (2.12)). Kinematically, sizeable values of 77 ~ 1 allow s 4 /m 2 ~ 1 that leads to significant 
non-factorizable corrections. In other words, collincar bremsstrahlung carries away a large part of initial energy. 
Since the (^-dependent and (^-independent Born level partonic cross sections have different energy behavior, the soft 
radiation has different impact on these quantities. 

Our analysis shows that two more factors are responsible for perturbative stability of the hadron level asymmetry 
First, one can see from Figs. 2 and 3 that both ^-dependent and ^-independent Born level cross sections take their 
maximum values practically at the same values of 77. Second, at fixed target energies, the gluon distribution function 
supports the contribution of the threshold region. In other words, sufficiently soft gluon distribution function makes 
the collinear gluon radiation effectively soft at the hadron level. In detail, the role of the gluon distribution function 
in perturbative stability of the azimuthal asymmetry in heavy quark photoproduction is discussed in Rcf. [14]. 

Another remarkable property of the azimuthal asymmetry closely related to fast perturbative convergence is its 
parametric stability. 3 Our analysis shows that the pQCD predictions for the cos 2<p asymmetry are less sensitive to 
standard uncertainties in the QCD input parameters than the corresponding ones for the production cross sections. 
For instance, changes of \xf in the range m c < pip < 2y / m 2 + Q 2 /4 affect the quantity A(Q 2 ) in charm production 
by less than 7% at p = 0.025 and £ < 4. For the (^-integrated cross section, such changes of piF lead to 30% 
variations in the same kinematics. We have also verified that all the NLO CTEQ5 versions of gluon density as well 
as the LO parametrization [31] lead to asymmetry predictions which coincide with each other with an accuracy of 
better than 1.5%. 

Parametric stability of the azimuthal asymmetry leads to the scaling: with a good accuracy the quantity A(p, x, £) 
in Eq. (1.6) is a function of three variables, so that 

A chaim (p, x, w A Bottom (p, x, (3.2) 

at the same values of p, x and £. To illustrate this property, in Fig. 9 we plot the asymmetry parameter A(Q 2 ) 
defined by Eqs. (3.1) as a function of £ at p = 10~ 3 and 10 -4 . Since the soft-gluon approximation is inapplicable to 
heavy flavor production at high energies, we give only the LO predictions for A(£). The LO results for the charm 



3 Of course, parametric stability of the fixed order results does not imply a fast convergence of the corresponding series. 
However, a fast convergent series must be parametrically stable. In particular, it must be p,R- and /iF-independent. 
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and bottom cases are plotted by solid and dash-dotted lines, respectively. One can see that both curves coincide 
with each other with an accuracy of better than 1%. 

It is also seen from Fig. 9 that, at the HERA energy (p ~ 10~ 3 and 1CP 4 for bottom and charm quark, respec- 
tively), the pQCD predictions for the AA in heavy quark leptoproduction are large and can be tested experimentally. 



B. Nonperturbative corrections 

Let us discuss how the pQCD predictions for azimuthal asymmetry are affected by nonperturbative contributions 
due to the intrinsic transverse motion of the gluon in the target. Because of the relatively low c-quark mass, these 
contributions are especially important in the description of the cross sections for charmed particle production [1]. 

To introduce kr degrees of freedom, k g ~ zp + kx, one extends the integral over the parton distribution function 
in Eq. (2.4) to fcr-space, 

dzg(z, n F ) -> dzd 2 k T f(k T )g(z, (i F ). (3.3) 
The transverse momentum distribution, f[kr), is usually taken to be a Gaussian: 

p~ k T / (k T ) 

^ = -^T < 3 - 4 > 

In practice, an analytic treatment of kr effects is usually used. According to [32], the fey-smeared differential cross 
section of the process (1.1) is a 2-dimensional convolution: 



dxdQ 2 dpq T dip 



/Q—k T /\K T ) do { 1 \ 

d2kT ^w, d X dQM P N QT d^ T -i kT )- (3 - 5) 



The factor \ in front of kr in the r.h.s. of Eq. (3.5) reflects the fact that the heavy quark carries away about one 
half of the initial energy in the reaction (1.1). 

Values of the fey-kick corrections to LO predictions for the cos 2tp asymmetry in the charm production arc shown 
in Figs. 4—6 and 9 by dashed curves. Calculating the fey-kick effects we use (fc T ) = 0.5 GcV 2 . At fixed target 
energies, fey-smearing for A(Q 2 ) and A{x) is about 20-25% in the region of low Q 2 and x, respectively, and decreases 
at large Q 2 and x. In the HERA range, expected values of the fey-corrections are systematically smaller (see Fig. 9). 

Analogous calculations for the case of bottom production are presented in Figs. 7 and 9. It is seen that fey-kick 
corrections to the 6-quark AA are practically negligible in the whole region of Q 2 and x. 



IV. CONCLUSION 



In this paper we have investigated the impact of soft-gluon radiation on the cos 2^3 asymmetry in heavy flavor 
leptoproduction. The NLL approximation provides a good description of the available exact NLO results for Q 2 < 
m 2 at energies of the fixed target experiments. Our calculations show that the azimuthal asymmetry is practically 
insensitive to soft-gluon corrections in this kinematics. We conclude that, unlike the (^-integrated cross sections, the 
cos2(/9 asymmetry in heavy quark leptoproduction is an observable quantitatively well defined in pQCD: it is stable 
both paramctrically and pcrturbativcly, and insensitive (in the case of bottom production) to nonperturbative 
contributions. This asymmetry is of leading twist and predicted to be about 15% at energies sufficiently above 
the production threshold for both charm and bottom quark. Measurements of the cos2<^ asymmetry in bottom 
production at HERA would provide an ideal test of pQCD. 

Data on the charm azimuthal distributions from the COMPASS and HERMES experiments would make it 
possible to clarify the role of subleading twist contributions. Our analysis shows that, in the low- a; region, the AA 
is sensitive to the gluon transverse motion in the target. At high x, the intrinsic charm contribution [33] to the 
structure functions may be significant [34,35]. In detail, the possibility of measuring the intrinsic charm content of 
the proton using the cos 2ip asymmetry will be considered in a forthcoming publication. 
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